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Abstract: Electromagnetic plane waves provide examples of time-dependent open string 
backgrounds free of a' corrections. The solvable case of open strings in a quadrupolar wave 
front, analogous to pp-waves for closed strings, is discussed. In light-cone gauge, it leads 
to non-conformal boundary conditions similar to those induced by tachyon condensates. A 
maximum electric gradient is found, at which macroscopic strings with vanishing tension 
are pair-produced - a non-relativistic analogue of the Born-Infeld critical electric field. 
Kinetic instabilities of quadrupolar electric fields are cured by standard atomic physics 
techniques, and do not interfere with the former dynamic instability. A new example of 
non-conformal open-closed duality is found. Propagation of open strings in time-dependent 
wave fronts is discussed. 
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1. Introduction 



While huge efforts have been devoted to finding phenomenologically viable static back- 
grounds for string theory, the study of time-dependent string backgrounds has only started 
to receive due attention recently. While a string implementation of inflation is still lack- 
ing, several examples of cosmological singularities have been studied [Q, |2|, ^ ^ ^, ^ , in 
first quantized string theory, raising strong doubts as to the validity of string perturbation 
techniques. Regardless of the issue of singularities, our understanding of string theory in 
time-dependent backgrounds is still very limited, and can only benefit from the study of 
solvable, if simplistic, cases. 

Gravitational plane waves are well-known examples of time-dependent backgrounds, 
free of string corrections 0, ^, They have received renewed interest recently in 

relation to the AdS/CFT correspondence, as they appear in a particular limit of AdS and 
other backgrounds They can be studied both at the first quantized level and in string 
field theory using light-cone techniques, even in the presence of Ramond backgrounds ||l2|. 
After fixing the light-cone gauge, the worldsheet theory acquires a mass gap proportional to 
the light-cone momentum , so that high energies on the worldsheet correspond as usual to 
long wavelength in target space. In fact, non-conformal two-dimensional sigma models can 
be turned into on-shell backgrounds by dressing them appropriately with dilaton gradients 
and fluxes in the light-cone directions ||l^ . Time-dependent gravitational plane- waves have 



been recently studied in flj, |15 



Less familiar perhaps is the fact that electromagnetic waves can be exact open string 
backgrounds, free of a' corrections. This was first recognized in the context of the Born- 
Infeld string where it was shown that a transverse displacement of a Dp-brane ^(xj), 
when supplemented by an electric field Fqi = is an exact open string background 
if ^ is an harmonic function of Xj. Upon thinking of <I> as the extra component of an 
higher dimensional gauge field Ap+i, this amounts to considering a configuration with null 
electromagnetic field strength F = di^dx^ A{dx^ — dx^), i.e. an electromagnetic plane wave. 
The case of a uniform, time-dependent or rather its T-dual version, was considered 
recently in 19 1. It includes the case of a monochromatic plane wave discussed long ago 



in [^]. In this paper, we consider the generalization of these exact string backgrounds to a 
time-dependent non-uniform plane wave accompanied by a constant and uniform magnetic 
field in the transverse direction, 

F = di^{x+,x') dx' A dx+ + ^Bij dx' A dx^ . (1.1) 
Here <I> is an harmonic function in transverse space, for the open string metric derived 



from S m, H, Gij = 6ij + {2'Ka'fBikBjk- In the hght -cone gauge, this breaks conformal 



invariance through the boundary condition 

d^X^ + p+d,^{X\ X+) + BijdrX^ = , a E {0, tt} (1.2) 

while preserving conformal invariance in the bulk, {d^ — d'^)X^ = 0. For general this 
implies a complicated non-linear and non-local coupling between the left and right movers. 
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This problem however becomes Hnear for $ quadratic in x*, corresponding to a quadrupolar 
electric potential. Our main goal in this paper will be to solve for the dynamics of first 
quantized open strings in this background. More generally, we shall consider open strings 
ending on two D-branes with possibly different quadrupolar electric potential and uniform 
magnetic fields. 

The non-conformal boundary deformation ( ^.10 ) is very reminiscent of deformations 



that arose in studies of tachyon condensation and background independent string field 



theory |2^, 24, 44, with two important caveats: (i) the deformation is taken on the 
(null) boundary of the strip, rather than the disk or the annulus and (ii) the "tachyon 
field" $ satisfies A<I> = in contrast to AT + 2T = 0. In particular, the potential $ is 
unbounded from above and from below and has no global minimum. This implies that, in 
the absence of a magnetic field, charged particles and open strings will be repelled away 
from the center of the quadrupolar potential. This is the direct analogue of the defocusing 
of geodesies by a purely gravitational plane wave. This kinetic instability does not mean 
that the background is unstable, on the contrary studies of the gravitational wave case 
have concluded to their stability pi 



Let us outline our main findings: 

i) For a critical value of the electric gradient, we find an instability towards producing 
macroscopic strings at no cost of energy: for that value, the tensive energy of an 
open string, which goes like its length square, is overwhelmed by the energy gained 
by stretching the string in the quadratic potential. The string then grows arbitrarily 
until it reaches the condensator that created the quadrupolar field in the first place 
and discharges it. This critical electric gradient is a non-relativistic analogue of 
the critical electric field strength in ordinary Born-Infeld electrodynamics. A simple 
elastic dipole model suffices to capture it. 

ii) In the absence of a magnetic field, the dynamical instability towards producing 
macroscopic strings is shadowed by the kinetic instability of charged particles in 
a quadrupolar electrostatic potential. One way to cure the latter is to switch on a 
constant magnetic field in the unstable directions, thereby confining the string into 
Larmor orbits. This is the familiar Penning trap configuration used in ion trapping 
studies. We show that the critical gradient instability is impervious to this change. 

iii) Another way to achieve confinement of charges in a quadrupolar electric field is 
to modulate the electric field at a resonant frequency with the proper modes of the 
electrostatic configuration: this is known as the Paul trap or rf (radiofrequency trap). 
We show that it manages equally well to confine charged open strings. We further 
compute the effect of a non-resonant adiabatic variation of the electric field. 

iv) While the one-loop open string amplitude is trivial due to the null nature of the 
electromagnetic field, this is no longer the case when the light-cone coordinate x"*" 
is compactified. Under this assumption, we compute the one- loop free energy for 
open bosonic strings, and find that it coincides with the overlap of boundary states 
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in the closed string channel. This gives a new example of open-closed duality in the 
non-conformal regime, beyond the ones found in |26, p7[| . 



v) We consider the propagation of open strings in a time-dependent quadrupolar electric 
potential. As usual for plane wave background there is no production of strings at zero 
coupling 1 29 1, however strings do get excited as they cross the electromagnetic wave. 



We compute the mode production in Born's approximation, and find as in |8| that 
impulsive wave profiles, unlike shock waves, transfer an infinite amount of energy into 
the string. We observe that the Bogolioubov matrix fails to incorporate an hidden 
degree of freedom of the background, for which the backreaction can be computed 
straightfowardly. We also discuss the adiabatic and sudden approximations. 

These results call for a number of further questions: • first and foremost, in analogy 



with the NCOS construction |3C] can one take a scaling limit towards the critical elec- 
tric gradient, and decouple the quasi-tensionless macroscopic open strings from the closed 
strings ? • Can a large limit of such D-brane configurations give rise to an holographic 



dual to gauge theories with time-dependent light-like non-commutativity |31, ^] ? • Can 
one take into account in a consistent way the condensation of these macroscopic strings ? 
More generally, can one take into account the backreaction of open strings with ^ 



onto the background (IT) self-consistently ? Is the backreaction finite or does it destroy 
the background irremediably ? Can one compute string production at finite string coupling 
? • Can one generalize this construction to other integrable boundary deformations, such 
as higher derivative Gaussian operators (f Xd^X, corresponding to excited string states, 
or Sine-Gordon-type boundary deformations ? • Can one understand further the relation 
between worldsheet RG flow and target space amplitudes ? What is the space-time in- 
terpretation of a boundary RG flow between two conformal fixed points ? • What are 
the implications of the critical electric gradient for the higher-derivative corrections to the 
Born-Infeld action ? • It has been observed [Q, |l8| that open strings in electric or null 
electric fields share many features with twisted closed strings at a spacelike or null orbifold 
singularity Is there a gravitational analogue to this critical gradient instability ? • 

Finally, at a far more futuristic level, one may wonder whether Bose-Einstein condensates 
of strings may be obtained using our string traps, or if the critical gradient instability of 
charged elastic dipoles may be one day be observed in QCD. We hope to come back to 
some of these questions in the future. 

The plan of the paper is as follows. In Section 2, we show that the configuration ( |1 . 1| ) 
gives an exact open string background, identify the solvable cases, discuss the motion of 
charged particles in this background, and how the transverse motion can be stabilized. In 
Section 3, we introduce an elastic dipole model and find a critical electric gradient at which 
it becomes tensionless. We show how adding a magnetic field or modulating the electric 
field stabilizes the kinetic instability, while preserving the phase transition. In Section 

4, we carry out the first quantization of open strings in a quadrupolar electric potential, 
possibly stabilized by a magnetic field. The one-loop amplitude is computed and open- 
closed duality is checked. Time dependent electromagnetic waves are considered in Section 

5, in the Born, adiabatic and sudden approximations. Our conclusions are contained in the 
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introduction. A mathematically rigorous derivation of the stability diagram can be found 
in Appendix A, a derivation of the modular property of the string partition functions is 
provided in Appendix B, and somewhat unwieldy normalization factors are gathered in 
Appendix C. 



2. Exact electromagnetic waves and ion traps 



Null electromagnetic fields F = Fi^dx'' Adx~^ are particularly simple as they provide exactly 
conformal deformations of the open string boundary conditions. Here we generalize the 
solution of [17| to include time dependence and a constant magnetic field^, and discuss how 
the light-cone gauge breaks conformal invariance on the boundary. We further describe the 
motion of charged particles in such backgrounds, and how their kinematical instability may 
be cured. 



2.1 Electromagnetic waves as exact open string backgrounds 

Let us consider the Abelian gauge field 

A = ^{x\x^)dx+ + BijX^dx^ , F = dA , 



(2.1) 



where ^{x^,x^) is an arbitrary function independent of x~, and Bij a constant antisym- 
metric matrix. One may show that A solves the equations of motion of Born-Infeld elec- 
trodynamics, 







where G 



fJ,U 



(detG,^)i/^G'^^Fp, 
r]^,y -\- {2T\-Q')'^F^fjrj^P Fyp is the open string metric |21, 22 



(2.2) 



ds^ = 2dx'^dx~ + Gijdx''dx^ + {2Tiaf \di^{x'^ ,x'')Y dx+ + d^^{x+ ,x'')Bijdx^ dx+ (2.3) 



under the only condition that ^ be an harmonic function in transverse space, 

1 



G 



di VG G'^dj]^{x\x+) = {) . 



(2.4) 



In this equation, the volume element does in fact drop since the transverse part Gij of the 
open string metric depends only on the constant magnetic field, Gij = 6ij + {2-Ka')'^ Bi^Bju- 
For two transverse directions, Gij is proportional to the identity matrix, so (0) reduces to 
d]'^ = 0, hence # has to be the real part of an holomorphic function. It is interesting to note 



that the open string metric (2^) is a gravitational plane wave in Brinkmann coordinates, 
together with an extra magnetic coupling. Note also that the energy momentum tensor 



associated to the configuration (LI) is given by 



(2.5) 



^The possibility to add a magnetic field while preserving exactness was observed in a discussion with A. 
Tseythn. 
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where G^" is the inverse of the open string metric appearing in ( ^.31 ) , and we subtracted the 
rest energy of the D-brane itself. The hght-cone momentum and energy 
therefore vanish irrespective of the profile Furthermore, one may check that the 

background ( |2.1| ) preserves half of the supersymmetries, namely those satisfying F+e = 0. 
In addition to solving the Born-Infeld equations of motion, one can argue that the 



background (|2.lD subject to condition ( |2.4|) is in fact an exact solution of open string 
theory to all orders in a'. Indeed, as the gauge field in ( |2.1[ ) is independent of x~, the 
embedding coordinate cannot be Wick contracted and remains classical. The sole 
effect of the magnetic term in the boundary deformation ^ A^dX^ is to renormalize the 
transverse metric from 5ij to Gij. The same reasoning as in |T^] then guarantees that there 
are no divergences coming from the electric term ^ <I>(x*, x'^)dx'^ as long as ^ is harmonic 
in transverse space. 



The Abelian configuration (2.1) may be generalized further to the non-Abelian U{N) 
case, by choosing a different electric field ^ for each generator of the gauge group: as only 
j4+ is non vanishing, the non-Abelian part of the field strength F = dA + [j4, A] remains 
null as it should. The magnetic field Bij should remain Abelian however. 

2.2 Light cone gauge and non-conformal boundary conditions 

In the presence of the § A^dX^ boundary deformation, the embedding coordinates of the 
string X*((T, r) are therefore still free fields on the worldsheet, 

{dl - dl)X'^ = , (2.6) 

but with the non standard boundary conditions (setting lira' = 1) 

d^X+ = (2.7a) 
daX' + di^a{X+,X')drX+ + BijOrX^ =0 at a = Ua , (2.7b) 

d^X- - ^^^a{X+,X') drX' = (2.7c) 

where we assumed different electric potentials, <l>o smd ^>i, at the two ends of the string, 
(To = and ai = tt, as in the case of a general non-Abelian configuration. By ( |2.6D X* can 
be decomposed in left and right movers: 

X\T,a) = f{T + a)+g\T-a) (2.8) 

The conditions ( |2.7D are compatible with choosing the light-cone gauge X+ 
upon which the condition ( 2.7b| ) becomes 



'a , 



d^X' + p+di^a{X+,X') + BijdrX^ =0 at a = aa 
corresponding to a boundary perturbation: 

J dr (^p+<l>{X+,X') + Is.^.X^-^) . (2.10) 
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As usual, choosing the hght-cone gauge breaks conformal invariance, but unhke gravita- 
tional waves, electromagnetic waves preserve conformal invariance in the bulk of the two- 
dimensional field theory. On the other hand, such non-conformal boundary deformations 
are reminiscent of the tachyonic perturbations considered in the context of background 
independent string field theory |^3|, with two important caveats: (i) since ^ is harmonic, it 
cannot be bounded from below or above, in contrast to the tachyonic perturbations usually 
considered; (ii) the deformation is taken on the boundary of the strip, and is inequivalent to 
a deformation <f dr$(x*) on the boundary of the disk since conformal invariance is broken. 
Nevertheless, it is interesting that electromagnetic plane waves can provide a real-time 
setting to study non-conformal boundary deformations. 

2.3 Solvable null electric fields and T-duality 



While the equation (2^) in general implies an inextricable non-local relation between the 
left and right-moving components (|2.8| ) of the field X*. we will be interested in cases where 
this relation becomes solvable. This is in particular the case when 

1) $a = {Aa)i{x~^)x'' is linear in x^, yet an arbitrary function of x"*": this corresponds 
to an electric field Fa = {Fa)i{x'^)dx^ dx~^ uniform in the transverse coordinates. The 
equation (|2.9| ) becomes a first order equation with source, 

d^X' + BijdrX^ +p+{Fa)i{x+) = at a = aa. (2.11) 

Conformal invariance is therefore unbroken, and the field X* can be computed easily 
by linear response. For i? = this case is T-dual to a configuration of two DO-branes 
following null trajectories Ya{x~^) where dY^/dx^ = —A\{x'^). Since this case has 
been considered recently in [|l9|, we will not discuss it any further. 



2) = {ha)ij{x~^)x^x^ /2 is quadratic in x*, and arbitrary function of x^: this cor- 
responds to an electric field Fa = {ha)ij{x~^)x^dx^ dx~^ which now has a constant 
gradient in the transverse coordinates x*. The equation ( |2.9D is still linear but in- 
cludes a mass contribution, in general time-dependent: 



d^X' + p+iK)^j{x+)X^ + -Bi.drX^ =0 at a = aa. (2.12) 



Conformal invariance is now genuinely broken. This will be the case of most interest 
in this paper, where we will show that the system exhibits an interesting phase 
transition at a critical value of the gradients {ha)ij- 

3) Clearly, the two cases above can occur simultaneously, for = {■Aa)i{x'^)x'^ + 
{ha)ij{x~^)x^x^ /2. This amounts to varying in time the location of the saddle point 
of the quadratic field, x*^^ = —{h~^y^Aj, and again could be computed by linear 
response from ( p. 12 ). 



It would be interesting to investigate the behaviour of open strings for null electric fields 
with a critical point of higher order, but this will not be addressed in this work. 
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Let us pause to discuss the T-dual version of the boundary deformation ( 2.12| ). Since 



the boson X are free fields in the bulk, they may be dualized into free bosons X*(r, a) = 
P{T+a) —g^{T — a). In terms of X\ the non-conformal boundary condition ( p.l2 ) becomes, 
after differentiating with respect to r (and assuming ha to be constant), 

dlX' + p+{ha)ijd,P = 0, (2.13) 

This corresponds to a boundary deformation J {h~^)ijdrX^drX^ dr/p^ by a massive opera- 
tor of the open string spectrum. Physically, it amounts to attaching a (possibly tachyonic) 
mass h^^ /p^ at the ends of the string. Notice that deformations by higher massive string 
states / m^fdrX^dJ^X^dT also give rise to integrable boundary conditions. This may be 
useful in evaluating the effective action for massive string excitations. On the other hand. 



we may consider another variation on our boundary problem (2.12), where r and a deriva- 
tives are exchanged: 

drX' + p+{ha)ijix+)X^ - ^Bijd^X^ =0 at a = aa. (2.14) 



Despite apparences, this boundary deformation is not T-dual to ( 2.12 ), although it is still 



integrable by the same techniques. Instead, for vanishing B one may integrate (|1|) with 



respect to r, and find that it describes open strings ending on infinitely boosted D-branes 
falling along the field lines of <I>. 

2.4 Charged particle in a null electric field 

As a warm-up, let us consider a relativistic particle of mass m and unit charge in an 
arbitrary null electric field A = ^{x~^ ,x'^)dx~^ satisfying the harmonicity constraint (p.^), 
together with a constant magnetic field. The equations of motion following from the action 

" 1 



S 



— (drX^'y -em' 



dr + Af.dX'' (2.15) 



read, after choosing the gauge e = 1, 

(dVdr2)X+ = 

{d^/dT'^)X' + di^ drX+ + BijdrX^ = (2.16) 
{d^/dT^)X- - di<^> drX' = 

The first equation can be integrated to X^{t) = Xq +p'^T. Taken together, the second 
and third equations imply that 

H = ^{Pif +p+p~ +p+<i>{X+,X')+m^ (2.17) 

is a constant, where = drX^ — BijX^ andp^ = dX^ /dr — ^ are the canonical momentum 
conjugate to X* and to X+. It is in fact the Hamiltonian constraint^ enforced by the 
Lagrange multiplier e in ( 2.15| ), hence the constant can be set to 0, thereby enforcing the 



mass shell condition. The motion thus becomes that of a non-relativistic particle in a time- 
dependent "electrostatic potential" $(x^,x*). The occurence of non-relativistic kinematics 
is a standard feature of light-cone qu antization. 

dimensionally correct Hamiltonian would be obtained by choosing e = 1/m instead, leading to 
H H/m. 
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2.5 Kinetic instability and trapping 

Now let us assume momentarily that no magnetic field is present. Since the potential $ is 
harmonic, there are no local minima, but only saddle points, hence the motion is expected 
to be unstable. This statement should actually be qualified: in two dimensions, a critical 
point of the potential will look as <I> = ^{z"'/n) with z = x + iy. Ordinary saddle points 
correspond to n = 2, <I> = ^(x^ — y^) and indeed the motion along y is that of an inverted 
harmonic oscillator. For higher order critical points however, the motion is easily integrated 
to z(r) = [n(n — 2)r]^/(^~") and therefore converges to z = for late times. We will not 
be able to solve this non-linear equation for the string for n > 2, hence we will not discuss 
it further. On the other hand, there are several known ways to cure the instability of a 
quadrupolar electrostatic potential: 

a) Penning trap: A well known scheme used in atomic physics to trap ions is to add 
a constant magnetical field to confine the unstable directions of the electrostatic 
potential. The simplest example involves 3 dimensions, with a cylindrically symmetric 



quadrupole potential V{x) 



2z'^), and uses a constant magnetic field 



along z to confine ions in the (x, y) plane. The equation of motion 
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(2.18) 



leads to the dispersion relation 

■2ujb + e){L0^ + 2ujb + e){2e 



- UJ 



(2.19) 



where the 3 terms correspond to the left and right polarizations in the (x, y) plane 
and the z direction, respectively. The motion is stable iff all roots are real, i.e. 6^ > e 
and e > 0. Hence only positively charged (say) particles are trapped, and only if the 
magnetic field is sufficiently strong. 

b) Paul trap: Charged particles can also be trapped with a quadrupolar field alone, 
if it is modulated at a particular frequency, $(x"'") = (e + /sin(a;r)(x^ — y^)/2, 
(parametrically) resonant with the proper mode of the quadrupolar field: this is 
known in atomic physics as the Paul, or rf (radiofrequency) trap. In mathematical 
terms, the equation describing motion along the x direction (say) is the Mathieu 
equation, 

x + [e + /sin(u;T)]x = (2.20) 



Its stability diagram, discussed in |37|, is qualitatively the same as for a rectangular 
modulation, which can be computed much more easily and is shown on Figure |3|. 
One sees that for e < 0, there is a small range of / where the motion is still stable. 
We will return to this point in Section 3.4. 

c) Quadrupolar trap: Finally, while a static quadrupolar electric field alone cannot trap 
neutral particles, it can trap polarizable molecules with negative polarizability: the 
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potential energy of an induced electric dipole //j = aEi is W = —aE"^; if a is 
negative, this is minimal at the center of the quadrupole, where Ei vanishes p8| . 
While the ground state of a molecule has usually positive polar izability, excited states 
if degenerate at E = can have negative a: the Stark energy shift due to the electric 
field is given by first order degenerate perturbation theory, and is negative for part 
of the multiplet. We will not pursue this case here. 

These schemes will prove convenient below to construct stable open string backgrounds. 
Note however that the instability encountered here is rather innocuous, as it is an instability 
of the motion of open strings rather than of the background itself. We will refer to it as a 
kinematic instability. Studies of a very similar issue in the context of gravitational plane 
waves with a non-positive quadratic form have found no hint of genuine instability [p5| . 



3. Elastic dipole in a quadrupolar electric field 

As we have just seen, relativistic particles of unit charge and light-cone momentum p"*" 
in the electromagnetic wave ( |l.lD behave as non-relativistic particles of charge in the 
(in general time-dependent) electrostatic potential ^ and uniform magnetic field Bij. In 
this section, we show that an analogous statement is true for a relativistic open string in 



(1.1). We consider a simple model of an elastic dipole in an electrostatic potential, which 



captures much of the physics of the string at low energies. 
3.1 Light-cone energetics and elastic dipole model 

Let us first call some basic features of light-cone quantization. As usual, the equations 
of motion (|2.6|) and (2.7) must be supplemented by the Virasoro conditions, which are 



impervious to the electromagnetic field. 



p-^d^X- + drX'd^X' = (3.1a) 
-p+drX- + i [{drX'f + {d^X'f] = (3.1b) 



Equation (3.1a) can be used to eliminate X in terms of the transverse coordinates, up to 



translational zero- modes Xq +Pq t. On the other hand, the Hamiltonian constraint ( |3.1b ) 
can be rewritten in terms of the canonical momenta 

7r_ = -drX+ = -p+ , m = drX' , 7r+ = -drX~ + $o(^(0))5o " (X(7r))<5^ , (3.2) 

and integrated from to vr to yield the mass-shell condition, 

p- =nic=^ [(vr.)' + id.X^f] da^ - + $i(X(7r)) (3.3) 



It is conserved if the potentials are independent of x"*". Comparing ( ^!^ ) to ( 2.17| ), we see 



that the light cone energy can be interpreted as the Hamiltonian of two non-relativistic par- 
ticles immersed in the time-dependent "electrostatic" potential j?"'"$i(X(7r)) —p"'"$o(^(0)), 
and bound by an elastic potential ^{do-X^)'^da. It will be crucial that the binding energy 
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grows like the length squared of the string, in contrast to a linear potential in the usual 
relativistic case. 

If we are interested with the low-energy aspects of the propagation of open strings in 
( |1.1| ), we may therefore content ourselves with a simplified model of a pair of charges, with 
positions x^(r) and x^(r), bound by a linearly stretched string, 

X\T,a) = ^[ax^j, + i7r-a)xl] (3.4) 

The light-cone energy of this "dipole" (it is a bona fide dipole only in the neutral case 
^o{x~^,x'^) = <I>i(x"^, X*)) is simply 

-P'^P~ + \{pI+Pr) + V{xl,xr) = (3.5) 

where V is the potential energy 

V{xl,xr) = ^{xi - x'jif+p+ [^i{xr) - $o(xl)] (3.6) 

Note that a static solution of this dipole problem will be an exact solution of the string, 
since ( |3.4D will then satisfy the bulk equation of motion {d^ — d'^)Xi = and the boundary 
conditions 

x''r-x'l - 7rp+ 9i$o(a;L) = at cr = (3.7a) 
x''r-x''l + vrp+ ai^>i(xj?) = at cr = vr (3.7b) 



In contrast to the single particle case, the two-body potential (3^) is no more harmonic, 
but satisfies 

(Al + Afi) y = — (3.8) 
vr 

where D is the number of transverse spatial dimensions. It may therefore admit stable 
extrema instead of only saddle points. 

3.2 Critical electric gradient and macroscopic strings 

We now consider the dynamics of an elastic dipole in a quadrupolar electrostatic field, 
= {ha)ijX^x^ /2. In the case where ho and hi commute, we can diagonalize them 
simultaneously and the potential energy ( |3.5| ) in a proper direction is then given by the 
quadratic form, 



i=k,l 



with in general an isolated degenerate point at the origin xl = xr = 0. It will be convenient 
to introduce the rescaled electric gradients Ca defined by: 

Ca = T^P^ha (3.10) 
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Figure 1: Left: Field lines of the quadmpolar electrostatic potential V = e{x'^ — 2/^)/2. Right: 
phase diagram for elastic dipoles in the (eo,ei) plane. The stable region is the domain in dark. 
Macroscopic strings get produced on the separatrix. 



Let us start by considering a single direction x. The stability of the potential V 
depends crucially on the sign of the determinant of Q, 

D = det{Q) = ei — eo — eoei (3-11) 

If (i) D > and ei — eo > —2, the two eigenvalues are positive and the origin is a global 
minimum. If (ii) D > and ei — eg < —2 the two eigenvalues are negative and the origin is a 
global maximum. Finally, in the intermediate region (iii) where D < 0, the two eigenvalues 
have opposite sign and the origin is a saddle point. This phase diagram is shown on Figure 
|l[ Note that in the limit of zero string tension, the only stable region would have been 
ei > 0, Co < 0, corresponding to a stable electrostatic well for the two ends of the string. 
Thanks to the elastic binding energy, the stability domain extends partly into the two 
quarters eoei > 0. If we now add a second direction y, due to the tracelessness constraint 
on h^, we are forced to choose opposite values (— eo, — ei) for the gradient along y. It is 
evident from Figure |l] that there is no way to have the motion stable along the two directions 
at the same time. More generally, due to the convexity of the domain (i), there is no way to 
choose gradients (eQ,e'^) along n directions in the stable domain (i) while still maintaining 
the tracelessness constraint Yli=i n(^0' ^1) ~ (0; 0)- This is consistent with the fact that the 
ground state of the string has a positive polarizability. One may ask whether choosing non- 
commuting electric gradients {ho, hi) at the two ends may somehow stabilize the system, 
however investigation of the crossed configuration $o = ^0(2;^ ~ ?/'^)/2,*I*i = hixy shows 
that this does not seem to be the case (see section 1.1.3| in the case of the open string). 



We will see in section 3.3 that stabilization can however be achieved by using a magnetic 
field, or by modulating the strength of the electric field. 

In spite of the fact that this electrostatic configuration is globally unstable as it stands, 
there is still a rather interesting phenomenon that takes place at the critical line D = 
0: there the potential admits a degenerate valley, corresponding to stretched strings of 
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arbitrary size L in the direction x, 



X = L{eoa — vr) , for ei — eo — eoei = (3-12) 

The physical origin of these states is clear: while the electrostatic potential decreases 
quadratically (when cqCi > 0) like the square of the length, the tensive energy of the 
non-relativistic string grows like the square length as well, and there exists a critical value 
of the electric gradient for which the two forces cancel for any size. When this value is 
reached, the strings become macroscopic and extend to infinity, thereby discharging the 
condensator that created the electric field in the first place. 

This phenomenon is very reminiscent of the critical electric field strength that appears 
in Born-Infeld electrodynamics, or equivalently for charged open strings in an electric field 



1 33, U]: the electric potential that pulls on the string ends grows with the length of the 
string L as EL, while the tensive energy grows as L/a' . At the critical value Ec = 
the two forces cancel and lead to strings of macroscopic size in the direction of the electric 
field, and small effective tension. If one attempts to go beyond Ec, the vacuum becomes 
unstable, tiny quantum strings get stretched to infinite distance very rapidly and discharge 
the condensator, thereby maintaining E < Ec- Such a regime has been used to construct 
the NCOS theory, a theory of non-commutative open strings decoupled from closed strings 



|3C]. It would be very interesting to see if such a decoupling limit could be taken here as 



well, leading possibly to a theory of non-relativistic open strings only (see |35] for other 
attempts to define non-relativistic limits of closed strings) . Another interesting question is 
whether and how these zero energy states can condense and backreact on the background 
(iLlp. This may be especially tractable due to recent progress in open string field theory. 

3.3 Elastic dipole in a Penning trap 

While the purely electrostatic configuration of the previous section exhibited an interesting 
phase transition, it still suffers of a global instability, that may shed doubt on the existence 
of these macroscopic strings. We now switch on a magnetic field, and investigate whether 
the same mechanism that allowed to trap charged particles still succeeds in stabilizing the 
string. We therefore consider our dipole model in the electrostatic and magnetic fields 

^a = ^ha(x^ + y^-Y^zA , B = bdxAdy (3.13) 



Here, we have taken into account the modified harmonicity condition (2^), which arises in 



dimension 3 and higher. According to our discussion in Section 2.1, this is the sole effect of 
a' corrections (we have set lira' = 1). Along the third direction z, the dynamics is exactly 
the same as in the purely electrostatic case up to a redefinition of (eo,ei), hence stability 
in this direction requires 

{l + b^ + 2eo) (1 + 6^ - 2ei) - (l + b'^f > (3.14a) 
(1 + 6^) +eo-ei >0 (3.14b) 
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In the transverse directions, the magnetic field leads to a modified dispersion relation, 



/l - eo - 

-2ibuj 1 



-1 



V 







2ibuj 


-1 



1 + ei - 





-1 

2ibuj 



VL 



-2ibuj 1 + ei 



7 



XvrJ 



[OJ] 



(3.15) 



At zero frequency, the effect of the magnetic field is nil, hence we still have production of 
macroscopic string states at the critical value of the electric gradients, 



D = ei — eo — eoei = 



(3.16) 



However, crossing this line does not lead to an instability anymore, as we now show. The 
discriminant of the dispersion relation (i.e. the resultant R of the polynomials d.et{Q) and 
d diet{Q) / duj'^) factorizes into 



[((l + 62)-eo) ((l + 62) + ei)-l] 



(3.17) 



When any of these factors vanishes, two roots uj^ of det((5) collide, and may, or may not, 
leave the real axis. Experiment shows that the third factor gives rise to an accidental 
degeneracy but no instability, while the fourth does indeed lead to the disappearance of 
two positive real roots into the complex plane, and hence to a genuine instability. The 
crossing of ( |3.16D only leads to the appearance of one zero eigenvalue which becomes 
positive immediatly again. If we restrict to the line eo + ei = for simplicity, we therefore 
have the following transitions: 

1. eo < 0: stable motion; 

2. eo = 0: one zero mode, macroscopic strings are created; 

3. < eo < h^: only real roots, stable motion; 

4. eo = 6^: two eigenvalues become imaginary; 

5. 6^ < Co < 2: partially unstable motion; 

6. eo = 2: one zero mode, macroscopic strings are created; 

7. 2 < eo < 2 + 6^: partially unstable motion; 

8. eo = 2 + 6^: two more eigenvalues become imaginary; 

9. eo > 2 + 6^: completely unstable motion; 

Here we have assumed that b^ < 2. If not, the stable region includes up to eo = 2, 
hence displays creation of macroscopic strings twice. In addition, there is an additional 
degeneracy at eo = 1 — 1/(46^), which can occur either in the stable or unstable region. 
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Figure 2: Phase diagram in the (eo,ei) plane for strings in a quadrupolar potential with small 
(left), intermediate (center) or large (right) magnetic field; in purple, boudary of stability domain; 
in red, line of production of macroscopic strings; in blue, line of accidental degeneracy. 



Examples of this phase diagram for several values of b are given in Figure |2|. To conclude, in 
the presence of a magnetic field, the stability region is enlarged to the upper left quadrant 
of the hyperbola: 



[l + b'^-eo) (l + 62 + ei)-l>0, 
eo < l + 



(3.18a) 
(3.18b) 



The stable region therefore includes all of the first critical line for creation of macroscopic 
strings, and extends into the formerly unstable region up to 



D 



2eo + Co ^2 



eo 



1 



for 6 < 1 



(3.19) 



where the kinematic instability takes over. If 6^ > 2, part of the second critical line for 
creation of macroscopic strings is included into the kinematically stable region as well. It 



is now possible to combine the two conditions ( |3.14| ) and (3.18) to find a region with global 
stability in all (x, y, z) directions. 

3.4 Elastic dipole in a Paul trap 



As we recalled in Sec 2.5, another standard way to trap ions in a quadrupolar electric field 
is to modulate the electric field periodically in time at a resonant frequency. Our aim in 
this section is to see to what extent the same mechanism can be used to trap charged 
elastic dipoles. Let us start by reviewing this mechanism in the point particle case. We 
thus consider a charged particle in a time-varying two-dimensional electric potential, 



y 



(3.20) 



where H{t) is a periodic function of t with period 2tt. After rescaling time, the equation 
of motion can be rewritten as 



X + [w^ + a^H{t)] X = 



(3.21) 
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where (a;^,a^) = ib(T/27r)^(/i, /) for X = x,y respectively. For H(t) = cost, this is known 
as the Mathieu equation, whose solutions are well studied transcendental functions. The 
stability diagram of ( p. 21 ) however depends little on the details of the function H{t), and 



can be studied much more easily in the case of a rectangular signal [p^], 

H{t) = 1 for t G [2n7r, (2n + l)7r[, , H{t) = -1 for t £ [{2n + 1)tt, (2n + 2)7r[ . (3.22) 
Since H(t) is periodic, the translation group must act by an element p of 5/(2, M) on a 



fundamental basis of solutions of the second order differential equation ( 3.21 ) , 




(t + 27r)= A\']{t), ad-bc=l, W{Xi,X2) = (3.23) 



[xl^ 

where W{Xi, X2) = XidtX2 — X2dtXi is the Wronskian, conserved in time. The eigen- 
values of this 5/(2, M) element can be written as e^^'^*'^ where the complex- valued Floquet 
exponent fi is computed from cos27r/i = (a + d)/2. The motion is stable iff fi is real, i.e. 
I cos(27r/i)| < 1. For the rectangular signal ( |3.22D , the matrix p can be computed straight- 
forwardly by matching X and its derivative at t = and t = tt. The Floquet exponent 
turns out to be determined by [^] 

1 f Xi X2\ 

cos(27r^) = cos2;icosx2 1 smxismx2 (3.24) 

2 \X2 XiJ 

where 



xi = vta/ lo'^ + ^ X2 = Tr\/ ij^ — , (3.25) 

The domain of stability | cos(27r^)| < 1 in the (a;^,a^) plane is depicted in Figure ^. For a 
small ripple a ^ u; and u;^ > 0, the motion is generically stable except when u; is integer or 
half- integer. At these values, the forcing is in resonance with the proper oscillation modes 
and an infinitesimal perturbation a is sufficient to destabilize the motion. This unstable 
region appears as a cusp on Figure |^, with boundaries well approximated at small a by 



1 



(-9+n2^ o 



= Y J^^, (3.26) 



2 f 2 ^,2 _L 7r^(2n+l)^-12 ^4 

2 ' ' ' 



1 -;K2^" + (2n+l)^^^ " J 



Conversely, an unstable motion at a = can be made stable by switching on a small 
perturbation. Indeed, there appears to be a stable region extending in the < domain 
above the line attaching at 0, well approximated by 

^2 ^ _!L«4 + J^^8 _ il!^!^«12 + 0(^16) (3 28) 

12 1512 9979200 ^ ' ^ ' 

It is therefore possible to choose the period and amplitude of the forcing (/, T) so that 
both (a;^,a^) and (— w^, — a^) are in the stable region, hence the motion is stable in both 
directions x and y. 
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Figure 3: Left: stability diagram for an harmonic oscillator with frequency modulated by a 
rectangular signal, in the (w^,q;^) plane. The stable region appears shaded, notice that it extends 
into the usually unstable regime liP' < 0. Right: zoom on the region where ^ 1, axes are 
exchanged. Instabilities set in at half-integer w, even with an infinitesimal amplitude of modulation. 



Let us now consider the case of an elastic dipole in a modulated quadrupole field. For 
simplicity, we assume that the same modulating frequency Tq = Ti = T is applied on either 
end of the open strings, with opposite amplitudes /o = — /i = /. The equation of motion 




1 — vrp^/iQ —1 
-1 l + irp+hi 









(3.29) 



can then be reduced to (3.21) for each of the proper directions of the constant matrix Q, 
with 



4 



a 



T 



27^p^ 

T 
27^p^ 



(3.30) 
(3.31) 



where the two signs correspond to the two modes of Q in the x direction. In the limit of 



small tension, the two eigenmodes uj^ become 

2 



UJ4 



27^p^ 



1 



-eo + l + 0(l//i)) , cu_ 



T 



27rp^ 



vr 



(ei + l + 0(l//i)) (3.32) 



The eigenmodes in the x and y directions are therefore not opposite anymore but shifted 
together upward. Nevertheless it is clear that for small enough tension one will still be able 
to choose the perturbation / such that both directions are stabilized. 

3.5 Elastic dipole in an adiabatically varying electric field 

After having discussed the motion of an elastic dipole in a resonant quadrupolar case, we 
now consider a dipole in a time-dependent quadrupolar potential ^a{x~^) = ha{x~^){x'^ — 
y^)/2, which we assume to vary adiabatically between two non-zero asymptotic values 
ha{±oo). In the particular case ho{x~^) = —hi{x~^) = h{x~^), the mass matrix Q{x~^) can 
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be diagonalized independently of time, yielding two decoupled harmonic oscillators with 
time-dependent frequency, 

i± + uj±{t)x± = (3.33) 

with uj^. = p^h{x^) lb 1, x-t = xl^xr. This problem can be solved for an arbitrary profile 
h{x~^) by separating the modulus and phase of, say, x± = s-i-e*'^±. The imaginary part of 
the equation relates the rate of phase variation to the modulus through Sj.'jj- = cte which 
can be set to 1. The real part gives a non-linear differential equation for the modulus, 

s±+w±(r)s± = 4- (3.34) 

The complete quantum mechanical S-matrix can then be found from the knowledge of s± 
at late times, where ( 3.34| ) reduces to the equation of motion deduced from the Hamilto- 



nian of conformal quantum mechanics |3^]. This result does not rely on any adiabaticity 
assumption, although determining s-|-(+oo) may not be feasible analytically. 

In the general case however, the proper directions of Q change with time, and it no 
longer helps to diagonalize the Hamiltonian. Instead, one may eliminate xl (say) and 
obtain a 4*'^-order differential equation for x = xr, 

TTX^*^) + (tr Q)x + det Q = (3.35) 

where we recall that tr Q = 2 + ei — eo and det Q = ei — cq — ei eo • It is not possible anymore 
to separate the phase and modulus of x = se^^ in general (although the imaginary part of 



(3.35) can be integrated once after multiplying by the integrating factor s(t)), however the 
problem can be solved in the adiabatic regime ha <ti h?^"^. Following the standard WKB 
method, we rescale t ^ tje and 7 — > 7/e. At leading order, ( |3.35| ) identifies the rate of 
phase variation with one of the instantaneous proper frequencies, 



7 = ^(r) = ±y ^ ± i^l + ^(eo + ei)2 , (3.36) 
hence 7(t) = uj{T')dT' , while the amplitude is given at next-to-leading order by 

{A{T)f/^ Jo T^Ar') 

where 

Air) = 2 (^(ei(r) - eoir) + 2)c^(t) - 2ioj{Tf^ (3.38) 

and uj denotes one of the branches of ( p.36| ) . The approximation remains valid as long as 
all eigenvalues stay separate from each other. In particular it breaks down at the line of 
production of macroscopic dipoles where two eigenvalues collide and leave the real axis. It 
would also break down if the electric gradient ha was switched off at ±00, except in the 
proper direction x^ which remains confined due to the string tension. 
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4. Open strings in a null quadrupolar electric field 



We now proceed with the first quantization of a string in a quadrupolar electrostatic 
potential = haijx'^x^ /2, where haij are constant traceless matrices, independent of 
x^. For now, we do not assume any particular commutation relation between /iq and hi. 

4.1 First quantization in a constant quadratic potential 

When decomposing X in left and right movers ( |2.8| ), the boundary condition ( |2.12| ) trans- 
lates into the system of linear equations 

f{r)-g{T) +p+ho{f{T)+g{T))+B{f{T)+g{T)) =0 
T^f(r)-g{T)+p+hi{T'fiT)+g{T))+B{T^f{T)+gir))=0 

where denotes the translation operator^, T'^f{T) = /(T-|-27r). In order to avoid clutter- 
ing, we have absorbed a factor p"*" into the gradients ha, and omitted the dxd indices. Such 
non-local differential problems have been studied in the mathematical literature under the 
name of linear differential difference equation of neutral type They are in contrast 

to the case of closed strings in gravitational waves, where the classical dynamics reduces 
to an ordinary differential equation for each mode of the Fourier expansion of the periodic 
coordinates X{a). They can nevertheless be solved rather straightforwardly using Fourier 
(or Laplace) analysis. 

4.1.1 Neutral string 

Before proceeding to the general case, let us first discuss the neutral case of an open string 
starting and ending on the same brane /iq = hi := h. Taking the difference of the two 
equations above, we get 

{T^-l){f+p+hf + Bf) = (4.2) 

Hence /(r) + p^hf{T) + Bf has to be a periodic function of r, with the Fourier series 
expansion (we assume [e,B] = 0), 

fir) + -/(r) + Bf{T) = ao + y({l + B) + —e) MnanC-'^^ = gir) - -gir) - Bg{T). 



IT ' ' ^ V nTT / vr 

(4.3) 

We chose a somewhat awkward normalisation of the a„'s in order to get the correct com- 
mutation relations at the end. ( f4.3D can be easily integrated to yield the normal mode 
expansion: 

e e 

-m{r+a) ^ m7r(l + ^) - 6 ^_^„(^„^)] (4-4) 

m7r(l — i?) -|- e 



where dxd matrices are still understood. The normalisation coefficients are given in 
Appendix The expansion therefore remains integer modded as in the e = case, 

^Its square root Tf{T) — /(r + tt) will be of use later on. 
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with the apparition of exponentially varying zero- modes (a+,a_). With the normalisation 
chosen, we have the standard commutation relations: 

K, a^] = n6'^n-m K, al] = i-e'^ (4.5) 

vr 

The Hamiltonian ( ^ ) can be easily computed: 

oo ^ 

p+W/c = a+a- + ^ 4a„ - — (4.6) 

n=l 

4.1.2 Charged string 

Let us now come back to the general case with cq ei, in the absence of a magnetic field 
for simplicity. Decomposing / and g into their Fourier modes, Equation ( |4.1| ) becomes a 
linear system 



IT TT 1 I V 



(4.7) 



The compatibility of this linear system puts a quantization condition on the energies w, 
which can be rewritten as the secular equation 

det [((vro;)^ + eico) sin(7rci;) — (ei — eo)vra; cos(7ru;)] = (4-8) 

Now let us assume that cq and ei commute (the case where they do not will be briefly 
considered in section 4.1. 3| ). The two matrices can then be diagonalized simultaneously 
and we can restrict ourselves to a potential of the form ^(x^ — y^) where x and y are two 
spatial directions. The problem thus decouples into two one-dimensional problems related 
by inverting (eo,ei) (— eo,— ei). The dispersion relation ([4. 81) thus reads 



2-Kiuj ^^'^ + eo «7ru; - ei 



or equivalently 



iTTOj — eo iiTUJ + ei 

tan(vr..) = / (4.9b) 

[TTLoy + eoei 

To find solutions of ( 4.9a| ) graphically, one has to look at the intersection points of the 
"Lorentzian curve" on the RHS with the tan curve on the lhs . The secular equation 
selects one real eigenmode u;„ in any interval [mr, (n + l)7r[, n E N*, with wave function 

X„(r, a) = Mn— ( e-^""(^+'^) + " go ^-.^„(.-.)\ ^ ^ ^ ^ ^^ ^^^^ 

where the normalization factor J\fn, given in Appendix ^ ensures that the commutators 
are indeed [an, aln] = uJnSn,m- 

The situation for the first segment [0, 1[ is more subtle however (see Figures ^ and |5|): 

• D = ei — eo — eoCi > and ei — eo > —2 

There is a single real root wo 7^ in ]0, 1[. As shown in Appendix]^, all the roots are 
then real and the motion is stable. The wave function Xq is the same as ( 4.10 ) with 
n = 0. 
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Figure 4: Graphical solution of the open string dispersion relation for a one-dimensional harmonic 
potential. Left: 13 > 0, ei — eo > —2, all roots are real. Right: _D < 0, ei — eo > —2, two roots 
have collided at to — and left into the complex plane. 




Figure 5: Graphical solution of the open string dispersion relation for a one-dimensional harmonic 
potential. Left: 13 < 0, ei — eo < —2, the root in the interval [0, 1] corresponds to the first excited 
level. Right: _D > 0, ei — eo < —2, two more roots have collided at w = and left into the complex 
plane. 



D = and ei — eo > —2 

As D approaches 0, the real root loq collides with its opposite —loq , leading to a pair 
of canonically conjugate modes with vanishing frequency. Their wave function can 
be obtained from the tOn ^ limit of ( |4.10| ) , and reads 

Xe(r, a) = M^il - ^a){xo + por) (4.11) 

TT 

where again the normalization (see Appendix^ gives the canonical brackets, [xcPo] = 
Those are exactly the macroscopic strings described in the previous section, occurring 
at a critical value of the gradient of the electric field. For small deviation away from 
the critical gradient, their energy is given by 



y 7r^(2 - 3eo + eg) 

D < (note that the line ei — eo = — 2 is contained in this domain) 

As D becomes negative, these two roots leave the real axis and become purely imag- 
inary, u)q = ±ikQ. This corresponds to unstable modes, with wave function 

X±(r, a) = f e±'=o(-+'^) + ~ ^° e^^°(— a± , [a+, a_] = iko (4.13) 

ko V =F7rfco + eo J 
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• D > and ei — cq < —2 

When crossing the D = curve in the region ei — eo < —2, the two roots ±wi cohide, 
hence two new macroscopic modes appear and then become unstable. They behave 
hke unstable modes described above in the case where D < 0. 

Altogether, we have thus found exactly the same phase diagram as in our dipole model, 
implying that the afore-mentioned instabilities have indeed a low-energy origin. Having 
found the normal modes, we can now carry out a standard mode expansion, focusing on 
the region of the stability diagram where ei — eo > — 2, 



n=l 



D>0 
D = 
D <0 



Xc 

x+ + x^ 



(4.14a) 



with canonical commutation relations. 



a+,a. 



(4.15) 



The contribution of the X coordinate to the light-cone Hamiltonian is easily evaluated in 
terms of these variables, and reads 



"Hx = ^ aian + < 



n=l 



D > : agOo 
D = 0: pI '^+Ex 

D <0 : a+a^ 



(4.16) 



The vacuum energy Ex has been evaluated using contour integral methods in [^]. The 
contribution of the boson Y can be obtained from ( [4.16| ) by simply reversing the sign of 
(eojCi) (note that D is not invariant under this operation). 

4.1.3 Non commuting electric gradients 



The dispersion relation (4.8) cannot be solved as easily as before, but we can examine the 
following crossed configuration: 



hn = h 



hi = h' 



(4.17) 



The secular equation still factorizes into 



■nuj 



V( 



2J2 



(4.18) 



The structure of roots is similar to the commuting case, except that the critical line is now 



ai D = + e"^ 



e^e''^ = 0, at which the slope at a; = of the lhs of ( 4.181) coincide the 



RHS (with the -|- sign). As before, for all n > there exists a real root in [n, (n + 1)[. For 
the branch with — sign, the slope at a; = will always be negative and there will never 
be an additional zero in the interval [0, 1[. However, there are always two complex roots 



- 22 - 



(Do = ±ikQ , that is to say, two unstable modes a+ and a_ (we use a " for oscillators of 
the minus branch, nothing for the plus branch). For the positive branch, the situation is 
comparable to the commutative case. If -D > 0, there also exists a real root ujq in [0, 1[. As 
D ^ 0, this roots collides with its opposite and leaves the real axis. For D < 0, it becomes 
a pair of complex conjugate roots loq = ±ikQ corresponding to the instability described in 
the previous section. At D = 0, we have macroscopic strings with arbitrary length, 



Yr 



3(e2 



2e3 



1-^ 



TT 



(xo + Pot) 



/e-1 



(4.19) 



and fixed angle Xc/Yc ■ 
configuration and reads 



Y^(e + l)/(e — 1). The Hamiltonian can be evaluated on this 



p+W = ^ alan + ^ alan 



+ < 



n=l 



n=l 



D > : agao 
D = 0: pI 

D <0 : a+a^ 



(4.20) 



Non-commuting electric gradients therefore lift part of the degeneracy of the commuting 
case, but do not prevent the kinematic instability to take place. 

4.2 Open string in a Penning trap 

While the sign of the quadrupolar potential was chosen so that the direction z was confining, 
the transverse (x, y) plane corresponds to unstable directions of the electric field. As we 
discussed in |3.3|, this can be stabilized by adding a constant magnetic field in the (x, y) 



plane. For simplicity we assume that e and B commute. The linear system (O) becomes 



— (1 + B^TTuj + eo 
(l + 5)«7ru; + ei)e-2^'^'^ 



^1 — B)muj + eo 
1 — B)mLV + ei 




The dispersion relation splits into two branches, 

7ra;(ei 



tan(7rw) 



eo) 



(1 + 62)(7ru;)2 ± 6(eo + ei)(7ra;) + eoei 



(4.21) 



(4.22) 



where the two signs correspond to the left and right-handed circular polarizations. Revers- 
ing iv exchanges the two branches, so that the two polarizations are canonically conjugate; 
we can therefore focus on the upper sign only. As before, a phase transition occurs when 
the tangents of the two curves on either side of (4.22) become equals. As b couples only 
at higher order in w, the transition line remains unaffected at D = ei — eo — eoei = 0. 
In contrast to (|4.9a ) however, only a single root collides with w = 0, as can be seen by 
expanding the dispersion relation at w = close to the transition line D — > 0, 



(eo - l)D 



(ttw) + 



(eo - 2)6 
eo 



:i-362)(3 



^^^±^{^ujf + 0{^^) = ^ (4.23) 



When Z) — > 0, one of the roots collides with 0, but simply crosses it as D changes sign. 
There is therefore no instability associated with the phase transition at -D = 0, instead 
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a creation of static macroscopic strings takes place as in the absence of a magnetic field. 
This is not surprising, since the effect of the magnetic field should only be felt in dynamical 
situations. On the other hand, an instability does take place as one goes further away from 
D = 0: indeed, by looking at the next term in the expansion (4.23), one sees that the two 
roots loq and uji collide and leave the real axis at a critical value 



4(eo-l)(3 



3eo + el] 



3e2(eo-2)2 



D 



D < 1 



(4.24) 



This is qualitatively the same as in our dipole model (3.19). As in the case of the dipole 



model in section |3.3| , the stability region is enlarged so as to include a part of the macro- 
scopic string creation line and to allow global stability region as regards x and y coordi- 
nates. Unfortunately, the precise critical line at finite D does not seem to be computable 
explicitely. At any rate, this demonstrates that the same set-up that allowed to stabilize 
elastic dipoles in a 3D quadrupolar potential also works for the full open string. 



4.3 One-loop amplitude and open-closed duality 

Just as gravitational plane waves exhibit no particle production nor vacuum polarization, 
the same is true for electromagnetic plane waves. A simple way to see this is to return to 
the Coleman formula for the one-loop vacuum free energy, written in light-cone variables, 

where we included the transverse momenta and string excitation levels in the label i, and 

denoted by the volume of light-cone directions. The integral over p~ enforces a delta 

function 6{tp~^), which implies that only p'^ = states contribute to the vacuum energy: 
but these states are insensitive to the presence of the null electromagnetic field, hence the 
vacuum energy reduces to that of flat space. Equivalently, in order to form a loop of open 
strings, one should make the worldsheet time r periodic, but this is incompatible with the 
light-cone gauge X'^ = Xq + p~^T unless p~^ = 0. 

On the other hand, it is well known that a non-trivial contribution to the vacuum 
energy is obtained when the light-cone coordinate is compact of radius R. Indeed, 



the integral over p~ in ( 4.25 ) is turned into a discrete sum on integers p~ = n/R, or after 
Poisson resummation into a periodic delta distribution with support at p^ = wR/t, w £ Z. 
One-loop vacuum energy with a loop modulus of 27rt writes 

^op=2^+- E X ^Trexp(ta2t)|^^^^^/^ (4.26) 



to=— oo 



The quantization of p"*^ is now in agreement with the fact that X^ should be a periodic 
function of r of period 27rt, modulo the identification X~^ ~ X~^ + 2itwR. For open strings. 
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= p^TLop and after a Wick rotation, one obtains 



^op = --F+_ / ^Trexp(-2^p+Hopt) 



\p+=Rw/t 

(4.27) 



to=— oo 



As usual, the open string vacuum amplitude should be reexpressible in the closed 
string channel as the propagator of closed strings between two boundary states describing 
the two D-branes with the electromagnetic flux on them. Since our D-branes have Neu- 
mann boundary condition along the light-cone coordinate , it is not possible to use the 
standard light-cone gauge to quantize the closed strings. Instead, one may use X+ = wRa, 
which follows directly from the open-string light-cone gauge = p^r after exchanging 
and rescaling the worldsheet coordinates, t ^ ta. The integer w is now interpreted as the 
winding number of the closed string around the compact coordinate . The closed string 
amplitude therefore reads 

A,i = Y,{Bi\ , _\ ^ \Bo) (4.28) 

where we used the fact that the boundary states Bq and Bi are annihilated by and . 
For closed strings propagation with Schwinger time nt 



oo 



oo 



oo 



Aci= yZ / di{Bi{w)\e-^'''^^\Bo{w)) 







(4.29) 

r-oo 

V / div+-z,i{t) 

^=-00-^0 

The equality between ( [4. 261 ) and ( |4.29D for D-branes in flat space now follows from the 
usual modular properties of open/closed string partition functions. 

We now return to the case of charged open strings in a constant quadrupolar elec- 
tromagnetic wave, = {ha)ijX^x^ and compute the one-loop vacuum energy in the case 
where the light-cone coordinate X+ is compact. Then, boundary conditions for open and 
closed strings involve parameter e and e which are 

Ca = T^P'^ha = TTwRha/t , Ca = TTwRha = TTp'^hat, (4.30) 

and thus obey the duality relation = eat (in fact, there is a i coming from the Wick 



rotation, see ( B.25 )). Starting with the open string channel, the contribution of the coor- 
dinate X with = cq and = ei on a cylinder of modulus q = e"^'^* is given, with due 
care paid to the zero-mode, by 



Zopit,eo,ei)=q^^l[{l-q' 



oo 

Ex TT n _ n'^'A^^ 



X < 

n=l 



(l_g^o)-i ifD>0 

HD^O (4.31) 



j_ / (2~3eo+eg) 
2t V 3(l-eo)D 
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In the zero field limit, one retrieves the correct behaviour 
Zop{t, eo, ei) ~ 



V 



\/eo - eit r]{t) ^/ki ri{t) Vt r]{t) 



(4.32) 



where V is the volume accessible to 0-modes and ri{t) = q^s ~ 9") the Dedekind 

eta function. The one-loop open string amplitude is obtained by putting together the 
contributions of the x and y directions as well as the remaining 22 transverse coordinates 
x*, and reads 



A 



op 



E 



w=—oo 



dt 



V22(47r^t??^(t)) ^^Zop{t, eo, ei)Zop{t, -cq, -ei) 



(4.33) 



where V22 is the volume of the 22 free transverse directions. For D < 0, the partition 
function picks an imaginary part, due to the unstable mode k^. It would be interesting to 
extract the production rate of macroscopic strings, in analogy with the stringy computation 
of Schwinger pair production in |40|. 

Now we turn to the closed string channel. As explained above, it is convenient to use a 
non-standard light-cone gauge X~^{T,a) = Rwa, compatible with the Neumann boundary 
condition on . Focusing on a single transverse boson X, the boundary condition ( 2.12| ) 
at (7 = becomes a condition at one cap r = of the cylinder. 



drX + -X = at r = 
vr 



(4.34) 



with e/vr given in ( 4.30| ). The coordinate X still satisfies the free field equation {d^—d^)X 
in the bulk, and can be expanded in the usual free closed string modes, 



Xci = xo + POT + i V - (a„e-^'^(^+'^) + a^e'^'^^^"'^; 



(4.35) 



The condition ( [4.34 ) can now be written as a condition on a boundary state \B{e)) in the 
closed string Hilbert space. 



1 + — On + I 1 

7rn J 



le 
Tin 



a- 



|i?(e)) = 



PQ + -XQ ) \B{e)) = 



(4.36a) 
(4.36b) 



which can be solved as usual by coherent state techniques, 
\B{e)) = M{e) e* 2e exp 

where 

M{e) = 



El nm + e 



\n=l 



n nrn — e 



|0,0) 



2\/2sinh(e) 



(4.37) 



(iH) 
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is fixed by open-closed duality as explained in Appendix |B|. Here we used a momentum 
representation to solve the zero-mode constraint (4.36b). The contribution of the coordi- 
nate X with Cq = eo and = ei to the one-loop amplitude in the closed string channel 
is now simply obtained as 

Z,i = J dpo{B{e,)\e-^'P'"^''\B{eo)) (4.38) 

where we included a summation over the closed string zero-mode po along the direction x. 
The free closed string Hamiltonian reads as usual 

1 °° 1 
^ ^ ^(„_^CK„ + a-nOn) - ^ (4.39) 

n=l 

The matrix element ( [4.38| ) can now be evaluated using the identity 

(O, 0| e^i^^e-^^^^^^+^^^^e^""^^^ |0, 0) = . (4.40) 

where we assumed the commutation relations [a, a^] = [q,q^] = 1. The integral on the 
zero-mode po is Gaussian. Altogether, we obtain, keeping in mind that the commutation 
relations for the a^s and the d^'s are [a„, at] = n and [d^, ali] = n : 

OQ / . I - • \ — 1 



\ 



t + i 



2 rr iim + eo nrn — ei 

U - -■ —-■ —Q 

1 L ) _~; V *7rn — eo mn + ei ) 



ei eo 



(4.41) 



where q = e . In the zero field limit, one obtains correctly 



Zciii eo, ei) ~ ~ J- ~ (4.42) 
Veo-eiT]{t) Vh7]{t) r]{t) 

One therefore obtains the one-loop amplitude in the closed string channel, 

°° r+oo 

Ad= Y] / dtV22^(t)"^^^c/(t,eo,ei)Za(t,-eo,-ei) (4.43) 
„ Jo 



w=--oo 



In order to establish the equality with the amplitude in the open string channel, we need to 
study the transformation law of Z under i — > A heuristic derivation can be performed 
as follows, generalizing techniques introduced in ||4^ for computing the vacuum energy - a 
full-fledged demonstration is provided in Appendix building instead on methods intro- 
duced in [27|. The logarithm of the open string partition function, disregarding vacuum 
energy contributions, can be written as the contour integral 

log Zop = ^ [ dz{log cl>e0^^^ (4.44) 

where dlog^o/ dz has single poles with unit residue at the open string eigenmodes, while 
log $c provides the correct contribution of that eigenmode to the free energy: 

^op{z) = 1 - e , ■ , ^d{z) = 1 - e ^^^^ (4.45) 

mz — eo iTTz + ei 
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The contour C follows the positive real axis from below from +00 — ie down to z = e < wq, 
and then back to +00 + ie from above. We may now unfold the contour to the imaginary 
axis C =] — zoo + e, —zoo + e[, and integrate by part: 

log = i- £ dz{\og 'J>o)^^^ (4.46) 

The term dXog^d dz now has single poles at z = in/t, corresponding to the closed string 
oscillators. The factor log $0 provides the correct contribution to the logarithm of the 
closed string amplitude, upon redefining cq = —iteQ,ei = —itei. More rigorously, we show 
in Appendix ^ the relation 

Zop(t,eo,ei)Zop(t, -eo,-ei) = Zci{i,eo,ei)Zci{i,-eo,-ei) (4.47) 

This gives a new example of open-closed duality in a non-conformal setting, beyond the 



computation in [26|. 



Finally, let us comment on the issue of fermionic strings in such backgrounds. The 
coupling of the Ramond-Neveu-Schwarz superstring to an electromagnetic background can 
be written as ^ dzip^^'^Ffj^i,. In the light-gauge, ip'^ = 0, hence the coupling vanishes 
for a null electromagnetic field This implies that the fermions remain quantized 

in integer or half-integer modes just as in flat space. While this seems to clash with 
worldsheet supersymmetry, the latter is broken by the choice of light-cone gauge^. This is 
another difference with tachyon condensation, where worldsheet supersymmetry requires 



to introduce a non-local interaction for the fermions [28|. The super symmetries unbroken 
by the electromagnetic background correspond to kinematical symmetries, which commute 
with the Hamiltonian and do not mix worldsheet bosons with fermions. 



5. Open strings in time-dependent electric field 

We now consider an open string moving in a time-dependent quadrupolar electric field. 
Due to the translational symmetry along x~ , there is as usual no production of strings at 
zero string coupling. On the other hand, a single string will generically get excited as it 
passes through the electromagnetic wave. This issue was studied in detail in for closed 
strings in gravitational waves. In this section we adapt their analysis to the case of open 
strings. 

5.1 Mode production and Bogolioubov transformation 

We assume a "sandwich wave" configuration, i.e. that the profile ha{x^) has a compact 
support in x^, and concentrate on a single direction x. At early times, we may expand the 
embedding coordinate X{T,a) into free open string modes, 

^ 2 

X(T,a) = xq + Pqt + i/ — o„ cos(n(T)e~*"'^ , r — > — 00 (5-1) 

^-^ n 

*We thank A. Tseytlin for emphasizing this to us 
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A similar mode expansion may be carried out at late times by replacing {xQ,po,an) by 
{yo,qo,bn)- Since the perturbation is linear, the two sets of modes will be related by a 
linear Bogolioubov transformation, 

f a (3 An 
7 S B„ 




n 



Bm B 







) 


Po 




\an) 



(5.2) 



Since the equation of motion ( [4.1D has real coefficients, this matrix satisfies the reality 
properties 

(a, /3, 7, <5) G MS {A, B, A, Bf^ = {A, B, i, 5)_„ , i3;;„ = (5.3) 

in addition to preserving the symplectic form —idxo A dpo + Yl,m>o'^'^rn A da-m/m. In 
particular, the off-diagonal components Bmn with m > and n < relate creation and 
annihilation operators at r ^ ±00, hence encode mode creation. This can be seen by 
evaluating the occupation number of the n-th mode at late times, 

I Z5 |2 



-"m,— n I 



where ip{x) is the quantum wave function of the open string zero-mode xq at early times. 

While the Bogolioubov transformation ( ^.21 ) suffices to describe the effect of the sand- 
wich wave on the string modes, it is important to observe that it fails to account for 
an extra degree of freedom which can be attributed to the background itself. Indeed, 
the zero-mode xq of the string is really the sum /o -|- go of zero-modes for the left and 
right-moving components, respectively. The difference /o — go usually has no physical sig- 
nificance, since it can be interpreted as a constant U{1) gauge potential A = {fo — go)dx 
along the direction x, which can be gauged away - indeed, under T-duality fo — go maps 
to the position of the D-brane along the transverse direction x, which is itself T-dual to 
the gauge potential A^. In a time-dependent situation however, the variation of fo — go 
between t = —00 and r = -|-oo has a gauge-invariant meaning, namely the integral of 
the electric field Ax{+oo) — Ax{—oo) = J F^^dx^ or equivalently the shift in transverse 
position of the T-dual D-brane. This change of the electromagnetic background induced 
by the propagation of strings is possibly the simplest instance of back-reaction in string 
theory. 

5.2 Born approximation 

Let us now consider the regime in which the electric field can be treated as a perturbation 
of the free open string case. This is in particular the case of an highly excited state of the 
string. We assume that the string starts out in a mode a„, and expand the left and right 
movers to first order in h, 

e"^"" + ( 7 1 (5-5) 
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where 6f and 6g are assumed to be of order h. The equation of motion ( |4.1| ) can now be 
rewritten as a Unear differential difference equation with source, 



D -D\ 6f\ ^ _2ian / eo 
L>r2 -d) [dgl nvr i (Tei) 



2 r.1 \ I - 1 M e"'"^ (5-6) 



where T is the shift operator by vr as above. The retarded Green function for this hnear 
system is easily computed by Fourier analysis 

and can be easily evaluated, with the result^ 

27r Vr + vrvr — r/ '^-^ \ 1 — 1 / 2'Kim 
\ / mez* \ / 

Note that G is the solution of the Green equation with a minus sign, VG = 5, with D, 
the matrix of differential operators on the lhs of ( |5.6| ) . One therefore obtains the solution 
after the electromagnetic wave has passed. 



-Y i /'"(eo-(rei))(p+r')^e-^("~'")^'dT'%-^'"^ 



vrn 



mez* 



(5.9) 



One deduces the matrix element of the Bogolioubov transformation 



• t-co 

B^n = 5mn + — (cq " (Tci ) ) (p+r)e-^("-"^)^dr (5.10a) 



oo 



= / (T(eo - (Tei))(p+r) - ^(rei)(p+r)) e-^^dr , (5.10b) 

Bn = — (eo - (rei))(p+T)e-*"^(ir , (5.10c) 

'— oo 



The other matrix elements can be obtained by starting with a zero mode xq + pqt at 



closely related computation was performed independently in [Q. 
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r ^ — oo , yielding 

1 r°° 

a = l-— (T(eo-(rei))(p+r)-7r(rei)(p+r))(ir , (5.11a) 



oo 



P = -^l {r\eo - {Tei))ip+T) - 7r{7r + 2T)iTei){p+T)) dT , (5.11b) 



7r2 



CO 

oo 



7 = A / (eo - (rei))(p+r)dr , (5.11c) 



oo 

oo 



1 

5 = 1 + —/ (T(eo - (Tei))(p+r) - 7r(rei)(p+r)) dr (5.11d) 



oo 



Am = ^ I (eo - {Tei)){p+T)e'^-dT , (5.11e) 



1 r°° 

^™ = / (^(«^o - (Tei))(p+r) - ^(Tei)(p+r)) e^^dr . (5.11f) 

J-oo 

and using the symplectic properties of the Bogolioubov transformation. One may in fact 
keep trap of the sphtting of the zero-mode 2:0 = /o + go between the left- and right-movers, 
and find the variation of /o and go separately between r = — 00 and r = +00: 

1 

Sfo - Sgo = — / eo{p+T)X{a = 0, r)dT (5.12) 

J-00 

As discussed in the previous section, this variation amounts to a correction to the null elec- 
tric background Fx+, and hence to a backreaction of the string on the electric background. 

As in the case of closed strings in time-dependent gravitational plane waves Q , it is of 
interest to ask if an open string can smoothly pass through a singular wave profile. This 
can be answered by computing the occupation number of the oscillator a„ after the wave 
has passed from (^.4|). If the wave profile /i" is a smooth (C°°) function with compact 
support in , its Fourier transform ha decreases faster than any power, and therefore the 
total excitation number is finite. In the case of an impulsive wave, where ha has a delta 
function singularity, the Fourier transform ha{u!) goes to a constant at \uj\ -^00: the total 
energy of the outgoing string ^ n{Nn) is infinite, implying that the string is torn apart 
as it goes through the singularity. For a shock wave with profile /i" proportional to the 
Heaviside step function, the total energy is instead finite, implying a smooth propagation. 
For the case of a "conformal singularity" h"'{x^) ~ l/xl discussed recently in the context 



of closed strings |15], the Fourier transform h"" diverges linearly and therefore the total 



excitation number itself is infinite. 



5.3 Adiabatic approximation 



We now consider the motion of open strings in an adiabatically varying electric field. The 
transverse coordinate X can still be decomposed into left-movers and right-movers. One 
may eliminate the latter to get an equation for /(r) only. 



= (1 - T^ )f 



'p^ %"'\ - T^) + -((Tei) - eo)(l + 

(rei)-eo TT 



/ (5.13) 
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+ 



^(eo - {Te,)T') - ^° (eo - {Te^)T') - -,eo{Te,){l - T^) 

TT vr (i ei j — eo tt 



where the dot denotes differentiation with respect to r, T is the shift operator by vr and in 
(Tei), T acts only on ei. The case eo = (Tei) =^ e gives a more simple equation 

(l_T2)/ + e(l-r2)/ = (5.14) 

which can be solved exaclty as an infinite sum of shifted solutions of y + ey = 0. However, 
the expression are hardly tractable and we will not use them. 

In line with the standard WKB approach, we assume that the phase of / varies much 
faster than its modulus. We thus write / = sexpi7, rescale 7 — > 7/e and time r — > r/e. 
Note that shift operator is also rescaled, Tx{t) = x{t + en). At leading order in e, we 
recover the standard dispersion relation 

{{T^iof + eoei) sin(7r7o) - (ei - eo)vr7o cos(7r7o) = (5.15) 



We choose for 7o(t) one of the branches in ( 5.15| ), which we assume to stay real at all 



times. This is for instance the case of the excited states Un>i- To next-to-leading order in 
e, we obtain the variation in amplitude (as well as a correction to the phase), 

so(r) = so(0) exp ( -ii, (7o(r) - 7o(0)) - f ^rr^.dA , (5.16) 



At') 



where 



A = — 7r7o - - [eo(l - eo) - ei{l + ei)] 70 + ^^(ei - eo - eoei), (5.17a) 

B = -27r7o7o + ip^ei^l + — (ei - eo)7o + 

vr 

+ - [eo(l - eo) - ei{l + ei)] 7070 + ^(eoei - eoei + ielei). (5.17b) 
vr vr 

This result is valid as long as e <C e^^^ and that no phase transition occurs. If a classically 
disallowed region is encountered, one should take into account the reflected wave using 
the standard WKB matching techniques. As in the dipole case, it does not seem possible 
to eliminate the phase altogether and get a deformed classical equation as in ( ^.34)) valid 
outside the adiabatic regime. 

5.4 Sudden approximation - open string in a Paul trap 

Finally, we consider the regime opposite to the adiabatic approximation, where the per- 
turbation is switched on and off by a Heaviside function of time. We have in mind the 
situation of the Paul trap, where the electric gradient is modulated by a square profile, 

$a = \hax^ , x+ e [nT, (n + ]^T)[ (5.18) 
^a = \h'ax\ x+ G [(n + ^T),(n + l)T[ . (5.19) 
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The embedding coordinate X can be expanded in modes a„ or as in (|4T4|) and (p^) m 
either of the two intervals. The two sets of modes can be related by matching X{a, r) and 
its derivative at = and = T/2 using the orthogonality relation ( |C.6D . At = 0, 
we find 



(5.20) 



where 



Snm{h' , h) = iuj'n / {X^n,h'dTXm,h — Xm,hdTX-n,h') da 

Jo 

= iuj'nJ^n,h'J^m,h 




with 01 



ITTUJn—ei 



(5.21) 



and the corresponding definition with primes on Cj for Q'^. When 



i7VLJn-\-ei 

h' = h, Snmih,h) = 5nm as it should. One therefore obtains the transition matrix from 
one period to the next, 



Tnm = ^ e'^'^^Snpih, h')e^'^'' Spmih' , h) 
pez 



(5.22) 



The issue of stability of open strings in the Paul trap is therefore reduced to the determining 
the modulus of the eigenvalues of the transfer matrix T. We expect a behaviour similar to 
the elastic dipole studied in |3.4|. 



Acknowledgments 

It is a pleasure to thank C. Bachas, G. D'Appollonio, M. Gutperle, B. Stefanski and A. 
Tseytlin for valuable discussions and suggestions, and especially C. Bachas and A. Tseytlin 
for critical comments on the manuscript. 



Historical note: After the first version of this paper was released on the arXiv, we 
became aware of |43, ^] where similar boundary conditions are discussed in the context 
of tachyon condensation at one-loop, with partial overlap with the mathematical results in 
Section 4. Section 4.3 on fermionic strings and 5.2 on neutral strings in time-dependent 
plane waves have been retracted in the present version, as it was realized that worldsheet 
supersymmetry need not be present on the light-cone, and that the appropriate condition 
Ci = T'^ei that leads to a simplification of the dynamics also renders the question of the 
Bogolioubov transformation moot. Finally, Section 4.1.3 on the one-loop string amplitude 
in the first release missed the fact, made apparent to us by the recent paper [^], that the 
light-cone needs to be compact in order to produce a non-trivial result; we have clarified 
this point in Section 4.3. As a reward to the reader for bearing with these errands, this 
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version provides added value in the form of an observation on an hidden degree of freedom 
of the background (Section 5.1 and 5.2), and a derivation of the open-closed duality formula 
(Appendix . 

Appendices 



A. Rigorous stability analysis 

In this Appendix, we establish the criteria for stability of open strings in a purely electro- 
magnetic quadrupolar trap on rigorous ground, following techniques explained in |4^]. The 
question to address is under which condition the characteristic equation A{z) = P[z,e'^) 
has all its zeros in the left half plane "^{z) < 0, where P{x, y) is a polynomial in x, y. 
For neutral differential difference equations, as is the case of interest in this paper, P{x, y) 
admits a principal term, i.e. contains a monomial arsX^'y^ such that all the other terms 
QmnX^y^ have either r > m, s > n, oi r = m, s > n, or r > m, s = n. The analysis is then 



based on the following theorem (see [41| for details): 



Theorem 1 Pontryagin (^^, A. 3) Let A(it/) = F{y) +iG{y) be the real and imaginary parts 
of A{iy) for j/ G R. If all zeros of A{z) have negative real parts, then the zeros of F{y) and 
G{y) are real, simple, alternate and 

G'{y)Fiy)^G(y)F'{y)>0 (A.l) 

for all 1/ £ R. Conversely, all zeros of A{z) will be in the left half plane provided that either 
of the following conditions is satisfied: 



1. (i) All the zeros of F{y) and G(y) are real, simple, alternate and (A.l) is satisfied for 
at least one t/ G R. 



2. (ii) All the zeros of F{y) axe real, and for each zero, (A.l) is satisfied 



3. (iii) All the zeros of G{y) are real, and for each zero, (A.l) is satisfied 



In order to establish condition (ii) or (iii), one rewrite or G as /(y, cos(y), sin(y)), which 
one decomposes as 

r s 

f{z,u,v) = Y,Y.'"'^'nl\n,v) (A.2) 

m=0 n=0 

where (jim^ is a homogenous polynomial of degree n in {u,v). The principal term in / 
is defined as the term z^4)i'^\u,v) such that either r > m, s > n, ov r = m, s > n, oi 



r > m, s = n in (A.2). We then have the theorem 



Theorem 2 Pontryagin (^^, ^.^^ If e is such that (fii''\cos{iy+e),sm{iy+e)) is non vanishing 
for all 2/ € R, then, for sufficiently large integers k the function F{z) = f {z , cos{z) , sm{z)) has 
exactly 4ks + r zeros in the strip ~2kn + e < K(z) < 2kn + e. Consequently, the function 
F{z) will have only real roots iff, for sufficiently large integers k, it has exactly 4fcs + r real 
zeros in the interval Ik = [—2kn + e, 2k-K + e]. 
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We now have the tools to investigate the stabihty of the hnear system Identifying 
z = 2Triuj, the characteristic equation can be written as 



A{z) = {-z^ + 4eoei)(e^ - 1) - 2(ei - eo)z{e' + 1) (A.3) 

Consequently, the real and imaginary parts of A{iy) read 

F{y) = {y^ + eoei)(cosy - 1) - 2(eo - ei)ysmy (A.4) 

G{y) = (y^ + 4eoei)siny + 2(eo-ei)y(l + cosy) (A.5) 



The condition (A.l) is satisfied trivially, since 

F'G -FG' = 2 + 4eoei) sin(y/2) + 2(eo - ei)y cos(y/2)] ' (A.6) 

as a consequence of the fact that the argument of A(iy) is y/2. The condition of stability 
is therefore that G{y) has only real roots. G{y) can be written as f{z, cos z, sin z) with / 
of the form ( |A.2| ), with principal monomial z^v. One may therefore choose e = '7r/2, and 
obtain the necessary and sufficient condition that G{y) should have exactly 4A; + 2 real 
zeros in any interval 1^ = [— 2/c7r + 7r/2,2A;7r + 7r/2]. The condition G{y) = consists of two 
branches, 

cos(y/2)=0 or tan (y/2) = ^^^^^ (A.7) 

y^ + 4eoei 

The first branch gives 2k zeros in the interval J^, leaving 2A; + 2 zeros to be found on 
the second branch. Here one meets the argument in the main text, and find that for 
D = ei — eo — eoci < 0, there are zeros for G{k) in 1^ {k large enough): this region is 
therefore unstable. For D > and ei — eo < —2, G{k) has only 4k — 2 zeros in Ik, showing 
also unstability. If -D > and ei — eo > —2 on the contrary, G has 4A; + 2 zeros in J^, hence 
has only real zeros by theorem 2. The motion is therefore stable. 

B. Open-closed duality 

In this section, we provide some details on the proof of the open/closed duality formula 
(4.47), taking inspiration from a similar discussion in [^]. Let's define Z by: 



^c,o(t, eo,ei) = Z^''°^'(t,eo,ei)Z^''°P(t,-eo,-ei) = Z^'°^(t, eo, ei)Z^''°P(t, eo, ei) (B.l) 

We start with the closed string one-loop partition function. Its logarithm can be 
written as 



logZ,(t,eo,ei) = (-27r£).2 ( j +2 log(AA(eo)AA(ei)) - ^ log ^ 



t + i 



1 1 

ei eo 



where q = e 



n=l 

'2lTt 




-. r^r — -g" + ^ -eo, ei ^ -ei . B.2 

tirn — eo iTrn + ei 
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Our aim is to perform a Poisson resummation on n. For this, it is convenient to add 
a regularizing mass x so as to let the sum on n runs from — oo to oo: 



n=l 



mn + eo iTrn — ei ^„ 
— Co «7rn + ei 



-— lim 

2 a;-*0+ 



i'n-Xn + ep - ei 
«7rx„ - eo ZTTXn + ei 

— log a; — log [277 it + i 



1 1 

ei eo 



where x„ = V?^^ + a;^. We now expand out the logarithm in the sum into 



(B.3) 



Si 



def 1 



iT^Xn + eo iTTXn - Cl 



n=— oo 

2 ^ L V^TTXn - eo iirxn + h ) ~^ K^Xn + eo «7rXn - ei 



^TTXn - eo ^7^Xn + ei 
^TTXn + eo ^7^Xn-el y _^ ^ ''^rXn - eo iT^Xn + ei 



^XnP 



p 



(B.4) 



The term q^"^/p term is best rewritten in integral form using a Schwinger time s, whereas 
the other factor can be expanded in power series with respect to ^ . The sum is now 



n,p,r,r' U 1 



(B.5) 



X^/ terms can be generated by differentiating the integral on s with respect to thus 
yielding an expression which can be Poisson resummed at ease: 



'1 O.A; ;;r' 



ds 



2V^ ^ , e?;'e?'^-'-' (dPy Jo s^^ \ w 

J_ 

27r ^ 



oo -j^ 

ds- 



^—p^s—vrs/t^—TT^Px^/s 



(B.7) 



The sum of h over Z can be folded into a sum over N at the expense of a factor of 2. Let 
us change variable s = expand e~^P^^ and differentiate r times with respect to P in S[ 
(the ' means: without n = term). We obtain 



1 roo °° 

n,p,r,r' l=U 



U T{1 + 
and for the n = term, 



Cq e^ 



27J- / ^ ~r' ,2r—r' 
p,r,r' 



dst 



{-Spy 1 



(B.9) 
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Gathering terms, one gets 

S[ = logZ, + 47riEo - log(47rx) - Tq - 2 log(AA(eo)AA(ei)) 



(B.IO) 



Let us now turn to the open channel. The logarithm of the open string partition 
function is 



log Zo{t, eo, ei) = -i27rt){Ex + Ey) log (l " ) " ^ V ~ ) ^^'^^^ 

p=0 p=0 

where q = e"^'^* and assuming that oJq'^ / 0. We denote by Px,y the set of solutions 
(renamed p) to the open string dispersion relation including w = solution. Introducing a 
regularizing mass x, the two sums can be rewritten (here we give only the expressions for 



J^log (l-q-v 
p=0 



lim^ ( 1 ~ 9 



p=0 



- lim V log (1 - q^i-) + - lim log (1 - o"^) . (B.12 



pePx 



2 x^o+ 



Note that we made a difference between the regularizing masses of closed and open strings 
because, as the duality involves a conformal factor j between the annulus and the cylinder, 



hence x = j. Therefore 



log(l-g'^) = log£ + log(27r) +0(2;) 



(B.13) 



We now expand the log and represent the power of g by a Schwinger-type integral, obtaining 
in the limit a; — s- 



1 



5/1/ -1 E .0.(1 -<,«)^^ E . 

We now represent the above sum by a contour integral 



ds 



1/2' 



(B.14) 



E 

pePx,PY 



dpe 



eo 



+ 



ei 



1 



iirp + eo «vr/9 — ei 

( 

ivrp — eo ivrp + ei 



1 + 



eo 



ei 



+ 



-2iiTp 



1 



TT^p"' + eg vr^'p"' + ef 

ivrp — eo ivrp + ei 

( 

ivrp + eo ivrp — ei 



-2iTTp 



(B.15) 



where C is a contour encircling all the p G Px,Py- C can be deformed into two lines Li 
and L2, where Li runs from 00 + ie to — 00 + ie above the real axis and L2 runs from 



—00 — ie to 00 — ie below the real axis. In order to expand the denominators of ( B.15 ) in 
power series of p^, one have to push Li and L2 to the region of the complex plane where 
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|/m(2;)| > eo, ei. The integral then picks contribution from the poles of the numerators of 
( |B.15| ), hence giving the following contribution 



2 \ / / • \ 2 



--log(l - e2^«)(l - e-2^°)(l - e^^'){l - e'^^^) 
2 ' (B.16) 



= — log (4sinh(eo) sinh(ei)) , 
where we anticipate the relation between and e^, (B.25). 



Expanding the denominators in ( [B.15| ), one can see that the contour integrals over Li 
and L2 are equal up to the p = term of the power serie expansion. This term gives the 
vacuum energy of the closed string computation: it reads 

;^l/:*f i?^-"--"'"'-'"'^ - ^ (4) - (-^'^) 

Terms with non vanishing p are on the other hand, summing Li and L2 equal terms, 



n, p6N* 



E / M 



iirp + eo ivrp ~ \ ^ ^ / ivrp — eo iirp + ei ^ ^ 
iirp — eo «vrp + ei / yivr/) + eo inp — ei 



eo 



ei 



7r2p2 _|_ g2 ^2p2 _|_ g2 



zvrp + eo 27rp — ei \ / inp — eo inp + ei 



ivT/o — eo ivrp + ei 



inp + eo i7rp — ei 



00 



s 



1/2' 



Changing variable to 



„2j0 ^ iTrp + eo iTrp - ei 
ivrp — eo iirp + ei 



(B.19) 



one finds the following expression, 
1 



^2 = — 1= / dp ( 2tt cos 2p(f> H — 

n,p 



i dcos2p(j)\ ^_2i^pp 



dp 



ds 







1/2- 



(B.20) 

where 5*2 is the sum without p = term. Integrating by part on p the second term, one 
gets 



s:. 



n,p 



ds 



1/2 



dp cos 2p<P ^e-"^----^*^(p^+-^)/^ 



(B.21) 



Let us now differentiate with respect to p, then expand 2cos{2p(j)) (with the same coeffi- 



cients as for the open string computation) and convert p^^ term to t^tf and the additional 
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p factor to ^ 

' in rt f v' " 



n,p,r,r 



Now we can integrate over p and differentiate with respect to p, 

n=l p=l r,r' U 1 ^ ' ^ ' 

(B.23) 

Expanding the exponential in power series to differentiate with respect to t^, one finds 



(B.24) 

and thus conditions of equality between 5( (closed side) and S'^ (open side) 

1 e- 

t = ^ , ei = = itci (B.25) 

There remains the vacuum energy term. We show that it is equal to Tq (substracting 
divergences) . 

J_ V — g-n2S-7r2t2(p2+^2)/5 ^ 1 >p }_ -2nht^/W+^ = 1 V - 27rtl»l^ 

sV2 2^n h^o 2 ^ \n '^7 

(B.26) 

where lim2;^o is assumed in the last line. On the other hand, the vacuum energy reads 

^ oo ^ 

p=o pePx,PY 
Thus vacuum energy is related to the term n = (substracting divergences) in the general 
computation. Then, performing the same manipulation than on the general h terms, and 
removing the linear and quadratic divergences arising in the definition of the vacuum 
energy, one can show that 

-{2Trt){Ex + Ey) = (term n = in the closed side), (B.28) 

Then 

S', ^ logZ„ + log ( ^^■■"■(^^o)smh(e.) ^ +'-fL-T, (B^29) 
From S[ = 82, one gets the following duality formula 

Zc{i, eo, ei) = {Af{eo)f{Af{ei)f (8sinh(eo) sinh(ei)) Zo{l,ieoi, iht) (B.30) 
Normalisation of the boundary states 

M{e) = ' (B.31) 

y'2V2sinh(e) 
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C. Normalization of eigenmodes 

C.l Without magnetic field 

• Normalization of the excited modes (n 7^ 0) and of mode when D > 0, 



47r 1 + 



ei 



eo 



(7ra;„)2 + ef (7rw„)2 + eg 



with 0: 



Normalization of the zero mode when D = 0, 



MS 



n{3{l-eo) + el) 
Normalization of the zero mode when D < 0, 



47r 1 + 



ei 



eo 



-(7rA;o)2 + ef -(7r/co)2 + e^ 



(Note that AT" = A/q replacing by iko.) 



(C.l) 



(C.2) 



(C.3) 



^ V^(l-e-2^^ 



(C.4) 



With the normalization chosen, we have the following orthogonality property between 
the Yn defined by 



Yn=Mn- 



iOr. 



ZTTLJn + eo 



(^^C^T^^ YmdrYfi) dcT — 5n,—n 
"^n 



(C.5) 
(C.6) 



C.2 With magnetic field 

The normalisation coefficient then become matrices M. We quote only the coefficients 
needed for the mode expansion ( |4.4| ). 

• Normalization of the excited modes (n 7^ 0) , 

1 



Mn 



Normalization of the zero mode for the two unstable modes: /cq = ij^^j 



(C.7) 



vr 1 - e i+s 



(C.8) 
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